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ABSTRACT 

Sgr A* is thought to be the radiative manifestation of a ~ 3.6 x 10 6 M Q supermassive black hole 
at the Galactic center. Its mm/sub-mm spectrum and its flare emission at IR and X-ray wavelengths 
may be produced within the inner ten Schwarzschild radii of a hot, magnetized Keplerian flow. The 
lightcurve produced in this region may exhibit quasi-periodic variability. We present ray-tracing 
simulations to determine the general-relativistically modulated X-ray luminosity expected from plasma 
coupled magnetically to the rest of the disk as it spirals inwards below the innermost stable circular 
orbit towards the "stress edge" in the case of a Schwarzschild metric. The resulting lightcurve exhibits 
a modulation similar to that observed during a recent X-ray flare from Sgr A*. 

Subject headings: accretion — black hole physics — Galaxy: center — magnetohydrodynamics — 
plasmas — Instabilities 



1. INTRODUCTION 

Sgr A*'s time-averaged spectrum is roughly a power 
law below 100 GHz, with a flux density S v oc v a . 
where a ~ 0.19-0.34. In the mm/sub-mm re- 
gion, howe ver, Sgr A* ' s spe ctrum is dominated by 
a "bump" (|Zvlka et alJ I1992D, indicative of two dif- 
feren t emission components ((Melia et al.l l2000t lAgoll 
l2000f) . High er frequenci e s correspond t o sm aller 
spatial scales (jMelial 11991 iNaravan et~all Il995h . so 
the mm/sub-mm radiation is likely produced near 
the blac k hole (BH). X-ray flares detected from 
Sgr A* dBaganoff et all [200lt [Goldwurm et al.1 12001 
iPorquet et al.l f2003: Belang er et al.ll2~005h may also have 
been produced within this compact region, either from 
a sudden increase in accretion accompanied by a reduc- 
tion in the anomalous viscosity, or from the quick ac- 
celeration of electrons near the BH (jLiu fc Melial [2002; 
iLiu et ail [2004). The energized electrons may also mani- 
fest themselves via enhanced emission in a hypothesi zed 
jet ([Markoff et al. 2001 . and references cited thereinl ). 

Near-IR flares detected from Sgr A* appe ar to be mod- 
ulated with a variable period sa 17 minute s (Gcnz el et al 
200a lEckart et aLll2006t iMever et al]l2006t lEckart et al 
2007f ). The X-ray and near-IR flares may be coupled via 



the same electron population, so one may expect simi- 
larities in their lightcurves. A long X-ray flare detected 
with XMM-iVewton in 2004 also appears to have a mod- 
ulated lightcurve, though n ot characterized by a fixed 
period ((Belanger et al.ll2008| ). If real, the modulation in 
both the near-IR and X-ray events is almost certainly 
quasi-periodic rather than periodic, with a decreasing 
cycle from start to end. But are the fluctuations due 
to a single azimuthal perturbation (i.e., a "hotspot"), 
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or from a global pattern of disturbance with a speed 
not dir ectly associated with the underlying Kep lerian 
period (|Tagger fe Melial l2006t iFalanga et al.1 12007( 1? In 
this Letter we examine the nature of the observerd quasi- 
period, and focus on its implications for the flow of mat- 
ter through the innermost stable circular orbit (ISCO). A 
principal result of this study is a ray-tracing simulation 
of the general-relativistically (GR) modulated lightcurve 
produced as the disrupted pl asma spirals inwards t o- 
wards the disk's "stress edge" (jKrolik fc Hawlevll2002( ). 

2. BACKGROUND 

Magnetohydrodynamic (MHD) simulations of Sgr A*'s 
disk demonstrate the growth of a Rossby-wave instabil- 
ity, enhancing the accretion rate for se veral hours, pos- 
sibly accounting for the observed flares (|Tagger fc Melial 
2006). The lightcurve produced by GR effects during 
a Rossby-wave induced spiral pattern in the disk fit 
the data relatively well, with a quasi-period associated 
with the pattern speed r ather than the Keplerian mo- 
tion ((Falanga et alJl2007t l. However, MHD simulations 
of black-hole accretion suggest that magnetic reconnec- 
tion might take place within the plunging region, due to 
the presence of a non-axisymmetric spiral density struc- 
ture, initially caused by the magnetorotational instability 
(MRI) associated with differen tial rotation of frozen-in 
plasma (see, e.g.. lHawlevll2001l ). 

In this case, the accreting flow is no longer Keplerian 
because of a radial velocity component. If Sgr A*'s quasi- 
period of ~ 17-25 minutes is associated with this kind of 
process rather than a pattern rotation, it would place the 
corresponding emission region at 0.73-0.94 nsco radii, 
below the ISCO (where risco = 3r s = QGM/c 2 ) for a 
Schwarzschild BH. Theoretically, we may therefore dis- 
tinguish the ISCO from the radius at which the inspiral- 
ing material actually detaches from th e rest of the mag- 
netiz ed disk — the so-called stress edge (jKrolik fc Hawlevl 
2002). The X-ray modulation would then be associated 
with the ever-shrinking period of the emitting plasma as 
it spirals inwards from the magnetic flare. 

Interest in "hotspots" began in the early 1980's in con- 
nection with quasi-periodic flux modulations observed in 
BHs accreting from a binary companion. The hotspots 
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are possibly overdense emission regions associated with 
magnetic instabilities. But even with a hotspot, a New- 
tonian disk does not produce a modulation since its 
aspect does not affect the total luminosity observed 
from it. Other than a dynamical periodicity (such as 
that due to an azimuthal, radial, or vertical oscillation), 
only GR effects can produce time-dependent photon tra- 
jectories resulting in a modu l ated lightcurve (see e.g . 
Cunningham fc Bardeenlll973t lAbramowicz et al.l 11991k 
Karas fc Bad 119921; IHollvwood et all 119951 : iFalcke et al l 
2000t iBromlev et all 1200 ll ). Even so, the "standard" 
disk picture of hotspot modulation has been based on 
Keplerian motion, for which one then expects a time 
variability directly related to the Keplerian frequency. 
Here, the modulation is not associated with such a fixed 
Keplerian frequency, but from a shrinking orbit and a 
monotonically decreasing period (see § ED . The relevance 
of hotspots has already appeared in (IHollvwood et all 
[19951 IMever et al.ll200a lEckart et al.1120071 : IMeliall2007L 
for review). What is lacking, however, is a non Keplerian 
treatment of the motion with the intent of probing the 
stress edge itself. 

So where exactly is the inner edge of the accretion 
disk in Sgr A*? This is a question asked in a broader 
context bv lKrolik fc Hawlevf (|2002D . whose MHD simu- 
lations of the plunging region in a pseudo- Newtonian po- 
tential identified several characteristic inner radii. Here, 
we assume a non-spinning BH, so our model pertains 
solely to the Schwarzschild case. 

The monotonic decrease of the period during the flares 
suggests that we are witnessing the evolution of an 
event moving inwards across the ISCO. The inflow time 
scale, ^inflow, which determines the rate at which plasma 
can move from one orbit to another, is given by t v = 
r g /v infiow w 9.6(r/r g ) 1 / 2 minutes (jliu fc: Melial 12002ft 
and is approximately 23.5 minutes at r — 3r s — 6r g , cor- 
responding to the ISCO for a non-rotating (i.e., a/r g = 0) 
BH. This time scale does not explicitly depend on a vis- 
cosity parameter since the viscosity is directly tied to 
the MRI physical p rocess via the induced Maxwell stress 
(|liu fc Me lia 2002). The inflow time scale defined here 
characterizes local processes occurring within the inner- 
most portion of the disk during the flares. By compar- 
ison, the dynamical time scale, td ~ 1-3 (r/r a ) 3 / 2 , i s 
roughly 19 minutes at this radius (|liu fc Melial [20021 
Thus, the azimuthal asymmetry giving rise to the mod- 
ulated flux during the flare may be due to a transient 
event associated with e ither a dynamical or viscous pro- 
cess close to the ISCO (jMelia et alll200lD . 

For a BH mass of 3.6 x 10 6 M Q , the inflow time scale 
at r ~ 2.5 r s (inferred from the average period) is just 
slightly larger than the average period, so the event could 
be due to the sudden reconfiguration of magnetic field 
lines frozen into plasma rapidly approaching the ISCO 
and then flowing across it towards the event horizon. 
Matter flowing past the ISCO may still remain "magnet- 
ically" coupled to the outer accretion flow, so a dynami- 
cally more meaningful radius is the so-called stress edge, 
where plunging matt er loses dynamical conta ct with the 
material farther out (Krolik & Hawley 2002T). This may 
simply be defined as the surface on which the inflow speed 
first exceeds the magnet osonic speed. 

In their simulations, Krolik & Ha wlevl (|2002p ) deter- 
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Fig. 1. — Upper panel: The stress edge radius, r s tressi in units 
of nsco, as a function of k, the exponent in the power-law formu- 
lation of Q(r). The dotted and dashed curves represent a period of 
17 an d 25 minutes, respec tively, using a black-hole mass of 3.6 X 10 6 
Mq JSchodcl et al. 2003]). Lower panel: The corresponding ratio 
of accreted specific angular momentum, j- ln , to the specific angular 
momentum at the ISCO. 

mined that this surface occurs somewhere between 
0.77nsco and nsco- The specific angular momen- 
tum j — r 2 il(r), in terms of the orbital angular fre- 
quency fi(r), continues to fall below nsco, though fl 
may not necessarily trace its Keplerian value, Quir) = 
(GM/r 3 ) 1 / 2 . In the absence of any magnetic coupling 
across nsco, matter would retain all of its specific an- 
gular momentum at the ISCO, so that the accreted 
value of j, which we will call j; n , would then simply be 
Jin = r isco ^if(nsco)- Instead, the MHD simulations 
show that j- m ss 0.95 j(nsco), for which r stress is then 
~ 0.8 nsco, within the range of values indicated by the 
location of the trans-magnetosonic surface. 

If the period in Sgr A* is decreasing monotonically, 
j(r) will not follow its Keplerian value below nsco- 
Therefore we will adopt the formulation f2(r) = flor~ K 
to fit the data in § 3. Clearly, k = 3/2 corresponds to Ke- 
plerian rotation; k is 2 in the extreme case of angular mo- 
mentum conservation. A reasonable fit to the data would 
therefore be associated with 3/2<k<2. At the bound- 
ary nsco, w e expect Q = £Ik, which then forces the 
constant £!o to have the value c^f^ r^ S( ?Q 2 . We calculate 
ntrcss using the quasi-periods 17 and 25 minutes emerg- 
ing from the X-ray lightcurve (see §[3|), and this is plotted 
as a function of k in Fig. [1] The radius r s t re ss falls within 
the range 0.73-0.96 nsco for all permitted values of k. 
The corresponding accreted specific angular momentum, 
for the same parameters as used before (see Fig. [TJ, is 
0.85 j(nsco) < jin < i(nsco) as a function of k. The ra- 
tio jin/j(nsco) = 0.95 from the MHD simulations would 
require k ~ 1.8, for which r s t r0S s ~ 0.77 nsco- These 
results are consistent with the MHD simulations, indi- 
cating that the infalling plasma below the ISCO remains 
magnetically coupled to the outer disk, though the dissi- 
pation of angular momentum is not quite strong enough 
in this region to force the gas into Keplerian rotation. 

3. THE INSPIRALING PLASMA MODEL 

With Q(r) known, we now incorporate strong gravita- 
tional effects in a geometrically and optically thin disk, 
describing the inspiraling disturbance using coordinates 
(r, 6, ip) in the co-rotating frame. In Fig. [2l we show 
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the inspiraling trajectory and duration of the emitting 
plasma. The observer is located at infinity with a viewing 
angle i relative to the z'-axis in the non-rotating frame, 
at (observer) polar coordinates (r', 9' , </?')■ The deflection 
angle of a photon emitted by plasma in the inspiraling re- 
gion is ip, varying periodically with cos ip = cos i cos <p, 
for a disk in the plane 9 = tt/2. Also, for G = c = 1, 
the BH's horizon occurs at r s — 2M, and the last stable 
orbit is located at rigco = 3r s . 




Fig. 2. — Upper panel: The inspiraling trajectory of the hotspot, 
beginning at nsco & n d terminating at 0.74risco- The dotted 
circle represents the location of the event horizon. Lower panel: 
The period as a function of the stress radius for the two extreme 
values of re adopted here, assuming a black-hole mass of 3.6 X 
10 6 M Q . 

We calcula te the lightcurve using a full ray-trac ing al- 
gorithm (see iLuminetl Il979t iFalanga et aTll2007t ). The 
disk from nsco to 90r s is an unperturbed, Keplerian 
flow, with angular velocity £Ik, and with specific angular 
momentum j'k = r 2 u v ju l — r 2 i7x. The corresponding 
four- velocity of the effective flow is then (it*, u r ,u ,u v ) = 
u* (1, 0,0, Cl K ), where u l = (1-MI /r)' 1 ' 2 (jMisner et all 
1973). The accretion flow is no longer Keplerian below 
the ISCO. 

Triggering a perturbation induces an azimuthal asym- 
metry in the region r strcS s ~ 0.73 < r < 0.9nsco- Below 
''ISCO; we use a simple representation of the bulk veloc- 
ity field, in which fl(r) = uf m /u{^, as described e.g., in 
iFukumura fc Tsurutal (|2004f) : 

< w = _A r e _(r " rBtrcBa:)/Asw sin 70 [fc I .(r-r s t rcS s)+"i( y 5/2-^ sw / 

. (1) 

In this case, the specific angular momentum is j m = 
r 2 uf w /ul w = £l r 2 ~ K . The subscript "sw" denotes 
the spiral wave, and the number m is the azimuthal 
wavenumber, fixed to be m = 1 for a single-armed spiral 
wave. The constant 70 = 2 is the width of the spiral 
wave, A r = 0.1 and A v = 0.1 are the amplitudes cho- 
sen to be relatively small, k r characterizes a tightness 
(i.e., the number of windings) of the spiral, and the ef- 
fective radial range of the spiral motion is controlled by 
A sw = 30, and ip sw = denotes the phase of the spi- 
ral. Since («g W , uf w ) is not axisymmetric, the net ve- 
locity field is also non-axisymmetric. For the effective 
flow then, (tt* w , u r sw , u e sw , uf w ) = u* w (l, < w , 0, fl r~ fe ), 
where u* w = [(1 - 2/m) - (1 - 2m/r)- 1 (< w ) 2 - r 2 ^] 1 / 2 , 
corresponding to the four- vector normalization condition 
g a ,pu a uP = -1. 

We consider four GR effects: (i) light-bending, (ii) 
gravitational Doppler effect defined as (1+z), taking into 
account the non-axisymmetric radial and azimuthal com- 
ponents below nsco : (hi) gravitational lensing, <iJ7 bs = 



bdbdip/D 2 (with D the distance to the source), expressed 
through the impact parameter, and (iv) the travel time 
delay. The relative time delay between photons arriv- 
ing at the observer from different parts of the disk are 
calculated from the geodesic equation. The first pho- 
ton arrives from phase ip = and r — nsco, and defines 
the reference time, To, which is set to zero. The observed 
time is then the orbital time plus the light-bending travel 
time delays, i.e., T ohs (ip 3vr ,r,i) = f2 _1 (7")</? sw + AT GR . 

The observed flux at energy E' is F ^ S (E') = 
Iobs(E')dQ b s , where I ^ S (E') is the radiation intensity 
observed at infinity and dtt b s is the solid angle on the 
observer's sky including relativistic effects. Using the 
relation I \, S (E' , a!) = (1 + z)~ 3 I enl (E, a), a Lorentz in- 
variant quantity that is constant along null geodesies in 
vacuum, the intensity of a light source integrated over 
its effective energy range is proportional to the fourth 
power of the redshift factor, / bs(a') = (l+z) _4 / m(^, <p), 
lem ( r if) being the intensity m easured in the rest frame 
of the inspiraling disturbance (jMisner etaj] [19731) . The 
disk radiates an inverse Compton spectrum, I cm , calcu- 
lated using the parameter scalings, rather th an their ab- 
solut e values. The spectrum parameters are (|Melia et al.1 
120011) the disk temperature, T(r), the electron number 
density, n e (r), the magnetic field, B(r), and the disk 
height H(r). This procedure gives correct amplitudes in 
the lightcurve, though not the absolute value of the flux 
per se. 

The synchrotron emissivity is therefore j s oc B n nt tx 
BT n ei where the nonthermal particle energy is roughly 
in equipartition with the thermal. The X-rays are pro- 
duced via inverse Compton scattering from the seed pho- 
ton number flux. Thus, with L soc d oc r 3 j s , where j s is 
the synchrotron emissivity in units of energy per unit 
volume per unit time, the soft photon flux scales as the 
emitted power divided by the characteristic area. That 
is, F sco d oc r 3 j s /r 2 — rj s , which is going to be roughly 
the same scaling as the seed photon density, so n scc d oc 
rj s oc rBTn e . The inverse Compton scattering emis- 
sivity is therefore ji C oc n nt n scc ^ oc (T n e ) 2 r B. Thus, 
jx ~ jio and the surface intensity is I cm oc J j x ds oc j x H, 
which gives finally I cm oc (Tn e ) 2 r B H. 
■ The flux at a given azimuthal angle if and radius r 
is calculated from a numerical computation of ip(a), fol- 
lowed by a calculation of the Doppler shift, lensing ef- 
fects, and the flux F bs as a function of the arrival time. 
For the persistent emission we use the best fit spec- 
tral parameters to t he Chandra data (jMelia et alj 120011 : 
iBaganoff et al.|[200lT ). described above as a surface emis- 
sivity I cm . The observed flare normalized flux is modeled 
with two polynomials, one between 0-100 minutes and 
the s econd from 100-160 minutes (see also I Meyer et al.1 
|2006() . The value k r is fixed at 11 to have the six observed 
cycles (see Fig. [3j solid line). The free parameters to fit 
the data are the inclination angle i and the k value. The 
integrated flux is calculated for an extended spiral wave 
90° long in the azimuthal direction and Ar = 0.28r ff in 
the radial direction, plus the persistent emission. The 
MHD simulations show that in the innermost part of the 
disk a spiral-arm often expands out to ~ 90° (see, e.g., 
Hawley 2001). The radial extent of the inspiraling region 
is set by the observed condition that six cycles should fit 
within the overall migration of the plasma from the ISCO 
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to the stress edge. In Fig. [3] (solid line), we show the best 
fit model for 72 ± 3° and k = 1.7 ± 0.05. 
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Fig. 3.— Light curve of the Augu st 31, 2004 flare in the 2-10 keV 
energy band (Belanger ct al. 2008), normalized with the observed 
mean count rate of 0.231 cts s _ l for the flare duration. The best 
fit model for an inspiraling disturbance is shown by the solid line 
using an inclination angle 72° and n = 1.7. The dotted curve 
represents a constant Keplerian period at the last stable orbit, i.e., 
r lSCO> * = 72°, and k = 1.5. Panels (a) and (b) show the residuals 
(in units of sigma) of the inspiraling and constant-period model, 
respectively, compared to the data. 



If we adopt the simple view that the last period cor- 
responds to the ISCO, then Sgr A* with a mass of 
3.6 x 10 6 Mq must be spinning at a rate a/r g > 0.2 — 0.4. 
With a more realistic analysis of the magnetic coupling 
between matter in the plunging region and that beyond 
the ISCO, we conclude that the peak of the instability 
probably occurs at ~ 0.97nsco> where the period is ~ 25 
minutes, and the flaring activity continues as the plasma 
spirals inwards, ending several orbits later when the mat- 
ter crosses the stress edge at ~ 0.8risco- 

The significance of the fit for an inspiraling distur- 
bance is x 2 /d-o.f. = 92.4/39, compared to \ 2 /d.o.f. = 
285.2/46 for a fixed Keplerian period (see dotted curve 
in Figure [3]). An inspiraling disturbance is preferred over 
a fixed orbit by a factor 2.6 in the reduced \ 2 ■ The resid- 
uals in the lower panels of Figure [3] show that the model 
using a fixed period produces modulations that are pro- 
gressively shifted in phase with respect to the data, by 
as much as ~ 16.5 minutes by the end of the flare. The 
inspiraling model, on the other hand, follows the evolu- 
tion of the flare and therefore fits the data much better. 
Plasma on such an orbit also produces a constant pulsed 
fraction= (I max -I m e a ,n)/(Im a ,x+Imin) of ~ 9%, compared 
with a linear increase from ~ 9% to ~ 11% for the in- 
spiralling wave; this effect is due to a radially-dependent 
gravitational lensing effect. Together, these two effects 
render the inspiraling scenario a better explanation for 
the data than the fixed orbit disturbance. 



4. CONCLUSION 
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